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ABSTRACT 

Sagittarius A* is the source of near infrared, X-ray, radio, and (sub)millimeter emission associated with the 
supermassive black hole at the Galactic Center. In the submillimeter regime, Sgr A* exhibits time-variable 
linear polarization on timescales corresponding to < 10 Schwarzschild radii of the presumed 4 x 10 6 Mq black 
hole. In previous work, we demonstrated the potential for total-intensity (sub)millimeter-wavelength very 
long baseline interferometry (VLBI) to detect time-variable - and periodic - source structure changes in the 
Sgr A* black hole system using nonimaging analyses. Here we extend this work to include full polarimetric 
VLBI observations. We simulate full-polarization (sub)millimeter VLBI data of Sgr A* using a hot-spot model 
that is embedded within an accretion disk, with emphasis on nonimaging polarimetric data products that are 
robust against calibration errors. Although the source-integrated linear polarization fraction in the models 
is typically only a few percent, the linear polarization fraction on small angular scales can be much higher, 
enabling the detection of changes in the polarimetric structure of Sgr A* on a wide variety of baselines. The 
shortest baselines track the source-integrated linear polarization fraction, while longer baselines are sensitive 
to polarization substructures that are beam-diluted by connected-element interferometry. The detection of 
periodic variability in source polarization should not be significantly affected even if instrumental polarization 
terms cannot be calibrated out. As more antennas are included in the (sub)mm-VLBI array, observations with 
full polarization will provide important new diagnostics to help disentangle intrinsic source polarization from 
Faraday rotation effects in the accretion and outflow region close to the black hole event horizon. 
Subject headings: black hole physics — Galaxy: center — techniques: inteferometric — submillimeter — 
polarization — accretion, accretion disks 



1. INTRODUCTION 

The Galactic Center source Sagittarius A* (Sgr A*) pro- 
vides the best case for high-resolution, detailed observations 
of the accretion and outflow region surrounding the event 
horizon of a black hole. There are several compelling rea- 
sons to observe Sgr A* with very long baseline interfer- 
ometry (VLBI) at (sub)millimeter 4 wav elengths. The spec- 
trum of Sgr A* peaks in the millimeter ( Markoff et al.ll2.Q07l 
and references therein). Interstellar scattering, which varies 
as the wavelength A 2 , becomes less than the fringe spacing 
of the longest baseline available to VLBI in the millimeter- 
wavelength regime. Indeed, VLBI on the longest baselines 
available at 345 GHz probes scales of twice the Schwarzschild 
radius (i?s) for a 4 x 10 6 M Q black hole. From previous ob- 
servations at 230 GHz, it is kn own that there are stru ctures 
on scales smaller than a few R$ dDoeleman et alj|2008l) . Such 
high angular resolution, presently unattainable by any other 
method (including facility instruments such as the Very Long 
Baseline Array), is necessary to match the expected spatial 
scales of the emitting plasma in the innermost regions sur- 
rounding the black hole and will be required to unambigu- 
ously determine the inflow/outflow morphology and permit 
tests of general relativity. 

This sensitivity to small spatial scales also makes millime- 
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ter polarimetric VLBI possible. Although the linear polar- 
ization fraction of Sgr A * integrated over the entire source is 
only a few percent (e.g.. [Marrone et ai1l2007l) . the fractional 
polarization on small angular scales is likely much larger. 
In general, relativistic accretion flow models predict that the 
electric vector polarization angle (EV PA) will vary along 
the circumference of the accreti on disk (Br omley et al.ll2001t 
iBroderick & Loebl 120051 120061) . indicating that single-dish 
observations and connected-element interferometers probably 
underestimate linear polarization fractions due to beam depo- 
larization. 

Polarized synchrotron radiatio n coming from Sgr A* was 
detected by lAitken et al.l ( 12000b at millimeter and submil- 
limeter wavelengths. Multiple observations since then have 
demonstrated that the polarized emiss ion is variable on 
timescales from hours to many days dBower et all 120051: 
iMacquart et"aT1l2006l IMarrone et al.ll2006al 120071 120081) . In 
one case, the timescale of variability and the trace of polariza- 
tion in the Stokes (Q,U) plane of the millimeter-wavelength 
emission are suggestive of the detection of an orbit of a po- 
larized blob of material (IMarrone et alj|2006bl) . Near infrared 
observations by iTrippe et al.l d2007h are also consistent with 
a hot spot origin for periodic variability. It is possible that 
connected-element interferometry may suffice to demonstrate 
polarization periodicity, but millimeter-wavelength VLBI, 
which effectively acts as a spatial filter on scales of a few to a 
few hundred Rs, can be more sensitive to changing polariza- 
tion structures. 

Initial millimeter VLBI observations of Sgr A* will nec- 
essarily uti lize non-imaging analys is techniques, for reasons 
outlined in lDoeleman et al.rd2009at henceforth Paper I). One 
way to do this is to analyze so-called interferometric "clo- 
sure quantities," which are relatively immune to calibration 
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errors dRogers et alJl!974l fl995l) . In Paper I, we considered 
prospects for detecting the periodicity signature of a hot spot 
orbiting the black hole in Sgr A* via closure quantities in 
total-intensity millimeter-wavelength VLBI. In the single po- 
larization case, it is necessary to construct closure quantities 
from at least three or four antennas in order to produce robust 
observables, since the timescales of atmospheric coherence 
and frequency standard stability do not permit standard nod- 
ding calibration techniques. Closure quantities can be used in 
full-polarization observations as well, but it is also possible to 
construct robust observables on a baseline of two antennas by 
taking visibility ratios between different correlation products. 
In this work, we extend our techniques to explore polarimet- 
ric signatures of a variable source structure in Sgr A*, with 
emphasis on ratios of baseline visibilities. 

2. MODELS AND METHODS 

We employ the same models discussed in Paper I to de- 
scribe the flaring emission of Sgr A*, and shall only briefly 
review these here, directing the reader to Paper I, and refer- 
ences therein, for more detail. These models consist of an or- 
biting hot spot, modelled by a gaussian over-density of power- 
law electrons, embedded in a radiatively inefficient accretion 
flow, containing both thermal and non-thermal electron popu- 
lations. 

The primary emission mechanism for both components 
is synchrotron. We model the emission from the thermal 
and nonthe r mal e lectro ns using the emissivitie s described in 
lYuan et al.1 (120031) and Uones & O'Delil (120071) . respectively, 
appropriately modified to acc ount for relativistic effects (see 
iBroderick & Blandfordl2004l for a more complete description 
of polarized general relativistic radiative transfer). Since we 
necessarily are performing the fully polarized radiative trans- 
fer, for the the rmal electrons we employ the p olarization frac- 
tion derived in lPetrosian & McTiernanl(11983l) . In doing so we 
have implicitly assumed that the emission due to the thermal 
electrons is locally isotropic, which, while generally not the 
case in the presence of ordered magnetic fields, is unlikely 
to modify our results significantly. For both electron popu- 
lations the absorption coefficients are determined directly via 
Kirchoff 's law. 

As described in Paper I, the assumed magnetic field geom- 
etry was toroidal, consistent with simulations and analytical 
expectations for magnetic fields in accretion d isks, though 
other field configurations are possible (e.g., iHuang et al.l 
120091) . While the overall flux of our models is relatively in- 
sensitive to the magnetic field geometry, the polarization is 
dependent on it. However, polarization light curves and maps 
with considerably different magnetic field geometries (e.g., 
poloidal) are qualitatively similar, showing large swings in 
polarization angle and patches of nearly uniform polarization 
in the images. 

Generally, synchrotron emission has both linearly and cir- 
cularly polarized components. However, the circular polar- 
ization fraction is suppressed by an additional factor of the 
electron Lorentz factor. For the electrons producing the mil- 
limeter emission, this corresponds to a reduction by a fac- 
tor of 30-100 in Stokes V in comparison to Stokes Q and 
U. Th is is consistent with observations by iMarrone et alj 
d2006al) . who obtain an upper limit of ~ 1% circular polariza- 
tion at 340 GHz. Therefore, we explicitly omitted the circular 
polarization terms in the computation of flaring polarization. 

In addition, we have neglected the potentially modest in- 
trinsic Faraday rotation. Within r < 10 2 -10 3 /?s the accreting 
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FIG. 1. — Integrated polarization traces of the models in the Stokes (Q,U) 
plane at 230 and 345 GHz over a full hot spot orbit, as would be seen by the 
SMA (for instance). 



electrons are expected to be substantially relativistic, and thus 
not contribute significantly to the rotation measure within the 
millimeter-emitting region. This is consistent with the lack of 
observed Faraday depolarization at the se wavelengths (e.g., 
lAitken et alll2000t IMarrone et al.ll2007l) . which itself implies 
the absence of significant in situ Faraday rotation. Similarly, 
beam depolarization caused by variations within an external 
Faraday screen on angular scales comparable to that of the 
emission region are empirically excluded. This leaves the pos- 
sibility of a smoothly varying external Faraday screen, which 
manifests itself in the VLBI data as an additional phase dif- 
ference between right and left circularly polarized visibilities, 
but does not affect our analysis otherwise. 

Model images are created in each of the Stokes parame- 
ters /, Q, and U. Six models differing in hot spot orbital 
period, black hole spin, and accretion disk inclination and 
major-axis orientation are produced at 230 and 345 GHz, as in 
Paper I. Model properties are summarized in TableQ] Source- 
integrated linear polarization fractions range from 0.8 to 26% 
for models including both a disk and a hot spot, depending on 
the model and hot spot orbital phase, with typical integrated 
quiescent polarization fractions (of the disk alone) of 10 to 
15%. Integrated EVPA variation over the course of the hot 
spot orbit ranges from 4 to 57°, depending on the model. The 
integrated polarization fractions and EVPA variations as well 
as the polarization traces in the Stokes (Q,U) plane (Fig.[T]i 
are all broadly consistent with the range of variability seen 
in the Submill i meter A rray (SMA) observations reported by 
IMarrone et alj d2006bl) . The local linear polarization fraction 
can be much higher, exceeding 70% in some parts of the ac- 
cretion disk. 

Simulated array data are produced by the Astronomical Im- 
age Processing System (AIPS) task UVCON for each of the 
Stokes parameters. The array is taken to consist of up to seven 
stations: the Caltech Submillimeter Observatory, James Clerk 
Maxwell Telescope, and six SMA telescopes phased together 
into a single station (Hawaii); the Arizona Radio Observa- 
tory Submillimeter Telescope (SMT); a phased array consist- 
ing of eight telescopes in the Combined Array for Research in 
Millimeter-wave Astronomy (CARMA); the Large Millime- 
ter Telescope (LMT); the 30 m Institut de radioastronomie 
millimetrique dish at Pico Veleta (PV); the Plateau de Bure 
Interferometer phased together as a single station (PdB); and 
a site in Chile, either a single 10 or 12 m class telescope (Chile 
1) or a phased array of 10 dishes of the Atacama Large Mil- 
limeter Array (Chile 10). Details of the method as well as 
assumed parameters of the telescopes are given in Paper I. 
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TABLE 1 
Model Parameters 
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Period 


i 


PA b 


V 


Disk Flux c 


Min Flux c 


Max Flux c 


Disk Pol. d 


Disk EVPA d 


Min Pol. d 


Max Pol. d 


Model 


(RgY 


(min) 


(°) 


(°) 


(GHz) 


(Jy) 


(Jy) 


(Jy) 


(%) 


(°) 


(%) 


(%) 


A 





27.0 


30 


90 


230 


3.19 


3.49 


4.05 


10 


-75 


7.0 


16 












345 


3.36 


3.63 


5.28 


1 1 


-81 


6.0 


24 


B 





27.0 


60 


90 


230 


3.03 


3.05 


4.03 


14 


-84 


6.9 


20 












345 


2.96 


2.99 


4.78 


13 


-80 


2.2 
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230 


3.03 


3.05 


4.03 


14 
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6.9 


20 












345 


2.96 


2.99 


4.78 


13 


10 


2.2 


26 
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0.9 


27.0 


60 


90 


230 


2.98 


2.99 


4.05 


15 


-86 


0.8 


21 












345 


2.96 


2.97 


4.00 


15 


-82 


1.6 


24 
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0.9 


8.1 


60 


90 


230 


2.98 


3.08 


4.15 


15 


-86 


10 


19 












345 


2.96 


3.04 


6.07 


15 


-82 


9.7 


24 


F 





166.9 


60 


90 


230 


3.07 


3.08 


3.38 


15 


-84 


9.8 


19 












345 


2.99 


3.00 


3.18 


13 


-80 


10 


17 



a Spin is given in units of the gravitational radius, Rq = GMc 2 = ■= R$ . 
b Accretion disk major axis position angle (east of north). 

c Stokes / flux density of integrated quiescent disk alone and minimum/maximum of system with orbiting hot spot. 

d Polarization fraction and EVPA of the disk emission alone and minimum/maximum polarization fraction of system with orbiting hot spot. 



3. POLARIMETRIC CONSIDERATIONS 

For ideal circularly-polarized feeds, the perfectly calibrated 
correlations are related to the complex Stokes visibilities 
(I u ,Q u ,U u ,V u )as follows: 

RR=1 V +V V 
LL = I V -V V 
RL=Q u + iU u 
LR = Q v -iU Vl 

where i = v— 1 and RL (for example) denotes the right circu- 
lar polarized signal at one station correlated against the left 
circular p olarized signal at another. We have used the con- 
vention of lCottonl (11993b . Other definitions, differing in sign 
or rotati on of the RL and LR ter ms by factors of i, are possi- 
ble ("e.g.. lThompson et al.ll200lh . but do not affect the analy- 
sis. Significant circular polarization is neither predicted in the 
hot spot models nor observed at the resolution of connec ted- 
element arrays dBower et alj|2003t iMarrone et al.ll2006ab . In 
the limit of no circular polarization (V„ = 0), /„ = RR = LL is a 
direct observable in the parallel-hand correlations, but Q v and 
U u appear only in combination in the cross-hand correlations. 
RL and LR visibilities, which are direct observables, are con- 
structed by appropriate complex addition of the Stokes Q v and 
U v visibilities. Right- and left-circular polarized (RCP and 
LCP) feeds are preferable to linearly-polarized feeds for de- 
tecting linear polarization, since the latter mix Stokes l v wi th 
Q v in the parallel-hand correlations (iThompson et al]l2001l) . 

For a point source, l v > \/Q^ + Ul + V^. However, for an 
extended distribution, the polarized Stokes visibilities can ex- 
ceed the amplitude of the Stokes /„ visibility. (For instance, a 
uniform total intensity distribution with constant linear polar- 
ization fraction but a changing linear polarization angle will 
produce no power in Stokes l v on scales small compared to 
the distribution, but the Stokes visibilities Q u and U v will be 
nonzero.) 

Analysis of polarimetric data is more complex than total 
intensity (Stokes /; we will henceforth drop the subscript 
on Stokes visibilities) data, but ratios of cross-hand (RL and 
LR) to parallel-hand (RR and LL) visibilities provide robust 
baseline-based observables immune to most errors arising 
from miscalibrated antenna complex gains. This stands in 
contrast to the single-polarization case in which robust ob- 



servables can only be constructed from closure quantities 
on three or more telescopes. The procedure for referencing 
cr oss-hand data to p arallel-hand data is explained in detail 
in ICottonl (1 1993b and lRoberts et al l (119941) and h as been used 



successfully in experiments (e.g.. IWardle 119711) . Several de 



tails warrant further discussion. We shall refer to the full ex- 
pressions for the observed correlation quantities: 

RR = R lR * = G lR G* 2R [ (In + Vn)e i( -^ 2) 

+D lR D* 2R (I l2 -V l2 )e i(+v '- ip2) 

+D ls P 2 V'' (+v,1+¥,2) 

+D* 2R P l2 e i( -^-^], 

LL = L X V 2 = G lL G* 2L [ (I 12 -V l2 )e i(+(p ^ 2) 

+D lL D* 2L (I l2 + V l2 )e i( -^ +ip2) 

+D lL P l2 e i( -^ 2) 

+D* 2L P 2 * x e i{+ ^ +V2 \ 

RL = R X L* 2 = G lR G* 2L [P 12 e i( - lp '-' P2) 

+D ls Z)? L P 2 V' (+tpl+V2) 
+D lR (In-V 12 )e i(+v "- ip2) 
+D* 2L (h 2 + V l2 )e i{ - lpi+V2) ] 1 

LR = L X R* 2 = G 1L G 2 * R [P 2 V' (+v,+lp2) 

+D lL D* 2R P l2 e i( -^-^ 
+D lL (In + Vn)e i{ -^ +V2) 

>2*( 



+D* 2R (In-Vn)e i(+v, - V2 \ 



where numeric subscripts refer to antenna number, letter sub- 
scripts refer to the polarization (RCP or LCP), a star denotes 
complex conjugation, G„x = g„xe'^" x is the complex gain in 
polarization^ £ {R,L} at antenna n, P = Q+iU, D„x is the in- 
strumental polarizatio n, and ip„ is the para llactic angle (equa- 
tions reproduced from lRoberts et al1ll994l) . The (p„ terms are 
constant for equatorial mount telescopes and can be incorpo- 
rated into the G and D terms, while for alt-azimuth mount 
telescopes the ip„ terms vary predictably based on source dec- 
lination, hour angle, and antenna latitude. It is likely that all 
of the telescopes in potential millimeter-wavelength VLBI ar- 
rays in the near future will have w,, terms varying with paral- 
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uv Distance (GX) 

FIG. 2. — Top: Visibility amplitude as a function of projected baseline 
length (v u 2 +v 2 ) for Model A at 230 GHz (noiseless). Stokes / is shown 
in black, and RL is shown in red. A real orbiting hot spot would persist for 
only a small fraction of a day, producing a plot corresponding to a subset of 
the above points. Contributions from the disk alone in the absence of a hot 
spot are shown in cyan (Stokes /) and green (RL). Bottom: Ratio of RL/I 
visibility amplitudes for the disk and hot spot (blue) and disk alone (orange). 
On small scales, RL/I can greatly exceed unity. 

lactic angle. 

The ratio of cross-hand to parallel-hand data (e.g., RL/LL) 
contains an additional phase contribution "J^, = V>nR _ V'nL 
equal to the phase difference of the complex gains of the 
right and left circular polarizations of antenna n dBrown et al.l 
1 19891) . These phase differences also enter into closure phases 
of cross-hand correlations as "fi + ^ + ^3. Fortunately, 
the right-left phase differences vary slowly with time (e.g., 
iRoberts et ai]|l994l) . since the atmospheric transmission is 
not significantly birefringent at millimeter wavelengths and 
both polarizations are usually tied to the same local oscil- 
lator. We will henceforth assume that the terms can be 
properly calibrated (for instance by observations of an un- 
polarized calibrator source), although proper calibration may 
not be strictly necessary for periodicity detection, since the 
expected timescale of variation of source structure is signifi- 
cantly faster than the timescale of variation of ^. Similarly, 
it is possible to determine the ratio of amplitudes of the real 
gains (r„ = g„R/g„L) from observations of a suitable calibra- 
tor. In general, r „ usually shows grea ter short-timescale vari- 
ability than \I/„ (IRoberts et al.l 11994b . Provided that proper 
instrumental polarization calibration is done, the fluctuation 
in r„ can be estimated from the RR/LL visibility ratio, since 
Sgr A* is expected to have no appreciable circular polariza- 
tion 5 . Even absent any complex gain calibration, it is prob- 
able that the contamination of the time series of cross-to- 
parallel amplitude ratios and (especially) phase differences 

5 Stokes V enters the expressions for RR and LL only as /± V, so even if 
circular polarization is detected, it will not prevent estimation of r„ unless the 
circular polarization fraction on angular scales accessible to VLBI is large or 
highly variable. 



by changes in r n and ^„ respectively will also be seen in the 
RR/LL amplitude ratio and RR-LL phase difference. Thus, 
large deviations seen in the cross-to-parallel quantities but 
not in the parallel-to-parallel quantities will likely be due to 
source structure differences, not gain miscalibration. 

Correcting for instrumental polarization (the D-terms) may 
be more difficult. Effe ctively, the D-terms m ix Stokes I into 
the RL and LR terms dThompson et al.l 120011) . Observations 
of calibrators with the Coordinated Millimeter VLBI Array 
(CMVA) at A = 3.5 mm found D-terms ranging from a few to 
21%, with typical values s lightly greater than 10% (lAttridgel 
120011: lAttridge et al] 120051) . Polarimetric observations with 
CARMA and the SMA in their normal capacity as connected- 
element interferometers have demonstrated that the instru- 
mental polarization terms on some of the telescopes that will 
be in cluded in future observations may be as low as a few per- 
cent (HowereLaDHQOl |MOTone!E3|2Q06i |2007l). How- 
ever, it is unknown how large the D-terms will be for poten- 
tial VLBI arrays at A = 1 .3 and 0.8 mm, as many of the critical 
pieces of hardware (including feeds, phased-array processors, 
and even the antennas themselves) do not yet exist for some 
of the elements of such arrays. In any case, contributions 
from the D-terms may be comparable to or larger than con- 
tributions from the source polarization, at least on the shorter 
baselines. The time scale of variations of D-terms is typically 
much longer than the time scale on which the source structure 
in Sgr A* changes, so carefully-designed observations may 
allow for the D-terms to be calibrated. At the angular res- 
olution of the SMA, polarization fracti ons of Sgr A* at 23 
and 345 GHz range between 4 and 10% dMarrone et al.l2007l) . 
although t he polarization fracti on may exceed this range dur- 
ing a flare dMarrone et al.l2008l) . Linear polarization fractions 
derived from single-dish and connected-element millimeter 
observations of Sgr A* are likely underestimates of the lin- 
ear polarization fractions that will be seen with VLBI, since 
partial depolarization from spatially separated orthogonal po- 
larization modes may occur when observed with insufficient 
angular resolution to separate them. That is, the small-scale 
structure that will be seen by VLBI is likely to have a larger 
polarization fraction than that observed so far with connected- 
element interferometery. 

Calibration of the electric vector polarization angle (EVPA) 
may be difficult, at least in initial observations, due to the 
lack of known milli meter- waveleng th polarization calibration 
sources (see, e.g.. lAttridgel l200"ll) . EVPA calibration will 
eventually be important for understanding the mechanism of 
linear polarization generation in Sgr A*, assuming that the 
linear polarization can be unambiguously corrected for Fara- 
day rotation. However, the ability of cross-hand correlation 
data to detect changes in the EVPA is unaffected by absolute 
EVPA calibration. 

In the low signal-to-noise (SNR) regime, the ratio of visibil- 
ity amplitudes can be a biased quantity. Visibility amplitudes 
are non-negative by definition, and the complex addition of 
a large noise vector to a small signal vector in the visibility 
plane will bias the visibility amplitude to higher values. Nev- 
ertheless, even biased visibility amplitudes may be of some 
utility in detecting changing polarization structure. Since the 
complex phase of noise is uniform random, phase differences 
are unbiased quantities. 

4. RESULTS 
4. 1 . Baseline Visibility Ratios 
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FIG. 3. — RL—I phase differences for selected baselines at 230 GHz. The same 2 hours of data, representing 4.5 periods (14.8 periods for Model E), are shown 
on all baselines except those involving PV or PdB. The solid line (red in the online edition) indicates the expected signal in the absence of noise, and the dots 
indicate simulated data at 8 Gbit s -1 in each polarization (16 Gbit s total). The blue line shows the EVPA that would be observed if the source were unresolved. 
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FIG. 4. — RL—I phase differences for selected baselines at 345 GHz. See FigurefJ]for details. 
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FIG. 5. — RL/1 amplitude ratio for selected baselines at 230 GHz. See Figure[3]for details. The blue line shows the integrated polarization fraction for the 
indicated model, shifted by 1 .0 for clarity. 




1.5 
1 

0.5 


1.5 
1 

.2 05 

K u 

° 1 s 

-d 1 ■ J 
=) 

£ 1 

Oh 

6 0.5 
< 

K 

K 1.5 

1 

0.5 


1.5 
1 

0.5 




PV-PdB 

Int. Pol. Frac. (+1) 



:: SMT — CARMA 





!'•'' ■ »' -s 



r 1 1 1 T f i ""pt 



0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 

Time (hr) 

FIG. 6. — RL/I amplitude ratio for selected baselines at 345 GHz. See Figures |5]and|5]for details. 
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FIG. 7. — Example polarimetric visibility ratio quantities with and without 
parallactic angle and D-term calibration. Noiseless model values are shown 
in red, simulated data (8 Gbit s per polarization) without parallactic angle 
and D-terms are shown in green, and simulated data with parallactic angle 
and D-terms are overplotted in black. Top: Simulated data for the CARMA- 
Chile 10 baseline of Model D at 230 GHz with zero D-terms. Parallactic 
angle terms effectively produce a small slope in the phase difference terms 
over timescales of interest and have no effect on visibility amplitudes in the 
absence of D-terms (i.e., the green and black dots are identical). Middle: The 
same baseline with 1 1 ± 3% D-terms. The inclusion of D-terms does not 
have a large effect because |P| is several times larger than |D/|. Bottom: The 
SMT-CARMA baseline with 11 ±3% D-terms. Uncalibrated D-terms will 
bias amplitude ratios and may produce qualitatively different phase differ- 
ences when the cross-to-parallel amplitude ratio is small but will not obscure 
periodicity even if the D-terms a re large (as is as sum ed in this worst-case s ce- 
nario based on the experience of Attridse 1 200 1 ) and Attridae et al. 1 2005) at 
86 GHz). The effect of uncalibrated D-terms will likely be largest on short 
baselines. 



While lower-resolution observations of Sgr A* find polar- 
ization of less than 10%, the effective fractional polarization 
on smaller scales can be much larger. Figure [2] shows the 
amplitudes of the (m,v) data 6 that would be produced by a 



disk and persistent, unchanging orbiting hot spot with param- 
eters as given in Model A at 230 GHz. The range in ampli- 
tudes reflects the changing flux density, both in total flux (i.e., 
the zero-spacing flux at u = v = 0) as well as on smaller spa- 
tial scales, as would be sampled via VLBI. Both total power 
(Stokes /) and polarization signatures fall off with baseline 
length, but on average the fractional polarization increases 
with longer baselines, and the ratio of Stokes visibility am- 
plitudes can exceed unity. All of our models produce much 
higher polarization fractions on small angular scales than at 
large angular scales, and all models except for Model F at 
345 GHz produce a substantial set of cross-to-parallel visi- 
bility amplitude ratios in excess of unity on angular scales of 
40-80 fias and smaller. 

We henceforth focus on ratios of cross-to-parallel baseline 
visibilities (e.g., RL/I). Plots of the RL-I phase difference 7 
are shown in Figures[3]and|4]for models at 230 and 345 GHz, 
respectively. At a total data rate of 16 Gbits" 1 , nearly all 
baselines exhibit signatures of changing polarization struc- 
ture. Due to the weak polarized signal on the longest base- 
lines, a phased array of a subset of ALMA (Chile 10, in the 
nomenclature of © may be required in order to confidently 
detect polarization changes on the long baselines, especially 
to Europe. The PV-PdB baseline (and to a lesser extent the 
SMT-CARMA baseline at 230 GHz) effectively tracks the ori- 
entation of the total linear polarization, since Sgr A* is nearly 
unresolved on this short baseline, and the calibrated RL phase 
of a polarized poin t source at ph ase center is twice the EVPA 
of the source (e.g.. lCottonlll995l) . 

Figures [5] and [6] show the RL/I visibility amplitude ratio 
for selected baselines at 230 and 345 GHz, respectively. The 
shortest baselines, PV-PdB and SMT-CARMA, effectively 
track the large-scale polarization fraction as would be mea- 
sured by the SMA, for instance. Because the short baselines 
resolve out several tens of percent of the total intensity emis- 
sion (as compared to the zero-spacing flux in Fig. |2]i but a 
much smaller fraction of the polarized emission, the variation 
in the RL/I and LR/I amplitude ratios is fractionally larger 
than in the large-scale polarization fraction. A bias can be 
seen in the amplitude ratios when the SNR is small, as noted 
in $3] (For brevity, we have shown only plots of the RL - 1 
phase difference and RL/I amplitude ratio. The LR-I phase 
difference and LR/I amplitude ratio exhibit similar behavior.) 

Closure phases of the cross-hand terms can be constructed 
in the same manner as for the parallel-hand terms, and these 
are robust observables. However, closure quantities are less 
necessary in the polarimetric case than for total-intensity ob- 
servations because robust baseline-based observables can be 
constructed. As Figure|2]shows, the visibility amplitude in the 
cross-hand correlations is much lower than that of the parallel- 
hand correlations on short baselines. The SNR of the closure 
phase is lower by a factor of y/3 than the three constituent 
baseline SNRs when the latter are all equal and is dominated 
by that of the weakest base line when there is a large differ- 
ence in the baseline SNRs dRogers et al.lll995l) . In Stokes /, 
the mean baseline SNRs (averaged over multiple orbits) are 
greater than or equal to 5 on virtually all baselines and all 
models at 16 Gbits" 1 total bit rate (8 Gbits" 1 each RR and 
LL) in a 10 s coherence interval, provided that the Chile 10 is 
used in lieu of Chile 1 . The number of triangles with SNRs 
greater than 5 on all baselines at 8 Gbit s" 1 in RL or LR is 



6 It is important not to confuse the antenna spacing parameters measured in 
wavelengths (conventionally denoted by lower-case u and v) with the Stokes 



parameters (denoted by upper-case U and V). 
7 Note that arg(RL) - arg(7) = wg(RL/I). 
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much smaller. Depending on the model, the SMT-CARMA- 
LMT and SMT/CARMA-LMT-Chile 10 triangles usually sat- 
isfy this condition, with Hawaii-SMT-CARMA also having 
sufficient SNR. Completion of the LMT, resulting in a system 
equivalent flux density of < 600 Jy at 230 GHz, will allow for 
strong detections on the Hawaii-LMT baseline and, impor- 
tantly, significantly strengthen detections on the LMT-Chile 
baseline. If the coherence time is significantly shorter than 
10 s, or if the obseved flare flux density is substantially lower 
than assumed in our models, closure phases may not have a 
large enough SNR to detect periodic changes. In any case, if 
polarimetric visibility ratios are successful in detecting peri- 
odicity, there may not be a need to appeal to closure quantities 
except insofar as they can be used to improve the array cali- 
bration. 

4.2. Instrumental Polarization Calibration 

We have also simulated the effects of not correcting for par- 
allactic angle terms and instrumental polarization by includ- 
ing Gaussian random D-t erms of 1 1 ± 3% wit h uniform ran- 
dom p hases, based on the lAttridgel d2001l) and lAttridge et alj 
( 12005b CMVA studies. This should be considered a worst- 
case scenario. D-terms at many of the telescopes will likely 
be subst antially better: e.g., 1-6% a t the SMA in observa- 
tions bv lMarrone etal.1 d2006al 12007ft and about 5% at the 6 
m antennas of the CARMA array dBower et al.ll20Q2b . These 
quantities affect the observed correlation quantities RR, LL, 
RL, and LR as indicated in §|3] We have ignored terms of or- 
der D 2 , but we have included terms of the form DP, since the 
polarized visibility amplitudes can be larger than the Stokes 
I=\ (RR+LL) visibility amplitudes on long baselines (Fig.[2j. 

Example data showing the effects of large uncalibrated D- 
terms is shown in Figure [7] Instrumental polarization adds 
a bias to the ratio of cross-hand to parallel-hand visibility 
amplitudes (e.g., RL/RR) as well as a phase slope and off- 
set to the difference of cross-hand and parallel-hand phases 
(e.g., RL-RR). These effects are much more pronounced on 
the short baselines, especially PV-PdB and SMT-CARMA, 
because the fractional source polarization on large scales is 
small (and thus \DI\ <jt. \P\). In most cases, the cross-to- 
parallel amplitude ratios and phase differences behave simi- 
larly whether instrumental polarization calibration is included 
or not simply by virtue of the fact that the polarized inten- 
sity is a large fraction of the total intensity. Deviations in 
the cross-to-parallel phase difference response appear qualita- 
tively large when the cross-hand amplitudes are near zero be- 
cause small offsets from the source visibility, represented as a 
vector in the complex plane, can produce large changes in the 
angle (i.e., phase) of the visibility. Large instrumental polar- 
ization can affect the expected baseline-based signatures but 
do not obscure periodicity, since source structure changes in 
Stokes / and P have the same period in our models. Of course, 
proper D-term calibration is a sine qua non for modelling the 
polarized source structure (but not for detecting periodicity). 
The D-terms can be measured by observing a bright unpolar- 
ized calibrator (or polarized, unresolved calibrator), and the 
visibilities should be corrected for instrumental polarization 
if possible. 

4.3. Periodicity Detection 

As in Paper I, we can define autocorrelation functions to 
test for periodicity. More opti mal methods exist t o extract the 
period of a time series of data (Rogers et al. 2009), but the au- 
tocorrelation function is conceptually simple and suffices for 



our models. The amplitude autocorrelation function evaluated 
at lag A; on a time series of n amplitude ratios A,- = RLi/U (or 
LRj/Ij) on a baseline is defined as 

j n—k 

ACF A (fc) ee V [(logAi-ziXlogAw*-/*)] , 

where /i and a 2 are the mean and variance of the logarithm of 
the amplitude ratios, respectively. The phase autocorrelation 
function is defined as 

^ n—k 

ACF^k) EE — - ^ cos(0 ; - <f> i+k ), 
1=1 

where </>, denotes the RL-I or LR-I phase difference of 
point By definition, ACF^O) = ACF (O) =1. The largest 
non-trivial peak corresponds to the period, with the caveat 
that the changing baseline geometries caused by Earth rota- 
tion can conspire to cause the autocorrelation function to be 
slightly greater at integer multiples of the true period. The 
phase autocorrelation function can suffer from lack of con- 
trast when the visibility phase difference is not highly vari- 
able as may be the case for the shortest baselines depending 
on the model (Fig. [8]), but the lack of contrast is not so severe 
as in the total-intensity case (cf. Paper I) due to the sensitiv- 
ity of short-baseline cross-to-parallel phase differences to the 
source-integrated EVPA. 

An array consisting of Hawaii, SMT, and CARMA is suf- 
ficient to confidently detect periodicity at a total bit rate of 
2 Gbit s" 1 (i.e., 0.25 GHz bandwidth per polarization) over 4.5 
orbits of data for Models A-E. This contrasts with the total in- 
tensity case, in which a substantially higher bit rate is required 
on the same array, depending on the model (Paper I). The key 
point is that the cross-to-parallel visibility ratios on the short 
baselines trace the overall source polarization fraction and 
EVPA, which are readily apparent even at the much coarser 
angular res olution afforded by co nnected-element interferom- 
etry (e.g.. iMarrone et ail l2006bl) . Long baselines are thus 
not strictly necessary to detect periodic polarimetric structural 
changes, although they will be important for modelling the 
small-scale polarimetric structure of the Sgr A*. In contrast, 
significantly higher bit rates and 4-element arrays are usually 
required to detect periodic source structure changes in total 
intensity (Paper I). 

Long-period models (e.g., Model F) are problematic for 
millimeter VLBI periodicity detection because it may not be 
possible to detect more than two full periods during the win- 
dow of mutual visibility between most of the telescopes in a 
potential VLBI array. As with the total intensity case (Pa- 
per I), the most promising approach for millimeter VLBI is 
to observe with an array of four or five telescopes, since 
large changes in cross-to-parallel phase differences and vis- 
ibility amplitude ratios tend to be episodic across most or all 
baselines. The LMT is usefully placed because it provides a 
long window of mutual visibility to Chile and the US tele- 
scopes as well as a small time overlap with PV, thus enabling 
a large continuous time range over which Sgr A* is observed. 
Assuming that the flare source can survive for several or- 
bital periods, connected-element interferometry may suffice 
to demonstrate periodicity. Sgr A* is above 10° elevation for 
approximately 9 hr from Hawaii and 12 hr from Chile. 

5. DISCUSSION 
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FIG. 8. — Autocorrelation functions for models A230 and C345 on selected baselines. Blue and red denote RL/I and LR/I, respectively. Green and magenta 
indicate RL/I and LR/I, respectively, using a completed 50 m LMT optimized for 230 GHz performance (fourth column) or Chile 1 (in place of Chile 10 in the 
fifth and sixth columns). The dashed line indicates the orbital period. 



5.1. Faraday Rotation and Depolarization 

VLBI polarimetry has several key advantages over single- 
dish and connected-element interferometry for understanding 
the polarization properties of Sgr A*. First, even the short- 
est baselines likely to be included in the array will filter out 
surrounding emission. Reliable single-dish extraction of po- 
larization information requires subtracting the contribution 
from the surrounding dust, which can dominate the total po- 
larized flux at 345 GHz and is significant even at 230 GHz 
dAitken et al.ll200Cil) . Contamination by surrounding emission 
is much less severe for SMA measurements, where the syn- 
thesized beamsize is on the order of an arcsecond, depending 
o n configuration (for in stance, 1"4 x 2. "2 in the observations 
of lMarrone et alJ2006al) . VLBI will do much better still, with 
the shortest baselines resolving out most of the emission on 
scales larger than ~ 1 mas (100 Rs), effectively restricting 
sensitivity to the inner accretion disk and/or outflow region. 

Second, the resolution provided by millimeter VLBI will 
greatly reduce depolarization due to blending of emis- 
sion from regions with different linear polarization direc- 
tions. Models of the accretion flow predict that lin- 
ear polarization position angles and Faraday rotation will 
be nonuniform througho u t the s ource dBromlev et al.l l200lt 
iBroderick & Loebl 120051 120061: iHuang et all 120081) . For 
this reason, ratios of cross-hand to parallel-hand visibilities 
(which are the visibility analogues of linear polarization frac- 
tions) can greatly exceed the total linear polarization fraction 
integrated over the source (cf. Figure [2] and Table [TJ. 

Third, VLBI polarimetry has the potential to identify 
whether changes in detected polarization are due to intrinsic 
source variability or changes in the rotation measure at larger 
distances. The former would be expected to be variable on 
relatively short timescales (minutes to tens of minutes), con- 
sistent with the orbital period of emission at a few gravita- 
tional radii. The latter would be expected to vary more slowly 
and affect only the polarized emission, not the total intensity. 
A cross-correlation between polarized and total-intensity data 



may allow the two effects to be disentangled. 

Linear polarization at millimeter wavelengths can be 
used to estimate the acc retion rate of Sgr A* (e.g., 
lOuataert & Gruzinovl 120001) . At frequencies below ~ 
100 GHz, no linear polarization is d etected due to Farad ay 
depolarization in the accretion region dBower et alj|1999allbl) . 
Linear po larization is detect ed toward Sgr A* at higher fre- 
quencies (lAifken et al. 20(3 e.g.,), where the effects of Fara- 
day rotation are smaller. Ultimately, accretion rates are con- 
strained by the lack of linear polarization at long wavelengths 
and its existence at short wavelengths. Measurements of 
the Faraday rotation exist, although it is unclear whether 
changes in detected polarization angles are due to changing 
source polarization structure or a variable rot ation measure 
dBower et alJl2005b iMarrone et al]l2006aLl2007l) . 

Longer term, imaging may be possible if all seven millime- 
ter telescope sites heretofore considered (and possibly others 
as well) are used together as a global VLBI array (e.g., an 
Event Horizon Telescope; iDoeleman et all l2009bl) . Imag- 
ing the quiescent polarization structure of Sgr A* may allow 
the characteristics of the source emission region to be distin- 
guished from those of the region producing Faraday rotation 
(which may overlap or be identical with the emission region). 
Contemporaneous millimeter VLBI observations at two dif- 
ferent frequencies would allow separate maps of the intrin- 
sic polarization structure and the rotation measure to be pro- 
duced. It may also then be possible to place strong constraints 
on the density (n e ) and magnetic fields (B) in Sgr A*. Briefly, 
the rotation measure is related to fn e B\\dl, while the total in- 
tensity is proportio nal to f n e B°! +1 d l, where a is the optically 
thin spectral index (Westfold 1959). Obtaining these results 
will require the ability to fully calibrate the data for instru- 
mental polarization terms and the absolute EVPA. It may also 
require higher i mage fidelity than a seve n-telescope VLBI ar- 
ray can provide (Fis h & Doelemanl2009l) . There are possibili- 
ties for extending a millimeter VLBI array beyond these seven 
sites by adding other existing (e.g., the South Pole Telescope) 
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or new telescopes (iDoeleman e t al. 2009bT). but full consid- 
eration of the scientific impact of potential future arrays is 
beyond the scope of this work. 

5.2. Physical Considerations 

The inclusion of a screen of constant Faraday rotation alters 
the phases of the cross-hand terms (and therefore the cross-to- 
parallel phase differences) but does not materially affect the 
detectability of changing polarization structure. The mean 
rotation measure of Sgr A* av eraged over multiple epochs 
is (-5.6 ± 0.7) x 10 5 radm" 2 (iMarrone et al.1 120071) . which 
corresponds to a rotation of polarization vectors by -55° at 
230 GHz and -24° at 345 GHz. Persistent gradients of rota- 
tion measure across the source are virtually indistinguishable 
from intrinsic polarization structure in the case of a steady- 
state source, but it is possible that the source structure and 
rotation measure change on different ti mescales, which woul d 
allow the two effects to be disentangled (IMarrone et al.l2007h . 
Comparison of changes in the polarization data with total 
intensity data (obtained from the parallel-hand correlations) 
and total polarization fraction (obtained from simultaneous 
connected-element interferometric data if available, else in- 
ferred from the shortest VLBI baselines) may be useful for 
identifying whether observed pola rization angle chan ges are 
due to a variable rotation measure dBower et al.ll2005l) . 

The physical mechanism that produces flares in total 
intensity and polarization changes is poorly understood. 
Connected-element interferometry at millimeter wavelengths 
has not been conclusive as to whether orbiting hot spots are 
the underlying mecha nis m that produces flares in Sgr A* (cf. 
IMarrone et alj|2006bl and IMarrone et al. 2008), or even as to 
whether multiple mechanisms may be responsible for flaring. 
Polarization variabili ty can be decorrelated from total inten- 
sity variability (e.g., iMarrone et aT]|2006bl) . and each shows 
variability on time scales ranging from tens of minutes to 
hours (and possibly longer). Spatial resolution will be key 
to deciphering the environment of Sgr A*, and thus there is a 
critical need for polarimetric millimeter-wavelength VLBI. 

Our results are generalizable to any mechanism producing 
changes in linear polarization, whether due to orbiting or spi- 
ralling hot spots, jets, disk instabilities, or any other mecha- 
nism in the inner disk of Sgr A*. Visibility ratios on baselines 
available for millimeter-wavelength VLBI will provide rea- 
sonably robust observables to detect changes in the polariza- 
tion structure on relevant scales from a few to a few hundred 
Rs, regardless of the cause of those changes. Clearly, period- 
icity can only be detected if the underlying mechanism that 
produces polarization changes is itself periodic, but baseline 
visibility ratios will be sensitive to any changes that are rapid 
compared to the rotation of the Earth. 

5.3. Observational Strategy 

Future millimeter- wavelength VLBI observations of Sgr A* 
should clearly be observed in dual-polarization unless not al- 
lowed by telescope limitations. Total-intensity analysis via 
closure quantities, as outlined in Paper I, can be performed 
regardless of whether the data are taken in single- or dual- 
polarization mode, but the cross-hand correlations can only 
be obtained from dual-polarization data. The cross-hand cor- 
relation data provide additional chances to detect variability 
via changing source polarization structure. 

By virtue of its size and location, which produces medium- 
length baselines to Chile and Hawaii as well as a long win- 



dow of mutual visibility with Chile, the LMT is a very use- 
ful telescope. The sensitivities assumed in this work for first 
light on the LMT may lead to biased amplitude ratios on the 
Hawaii-LMT baseline (Fig. |5), but this will not prevent de- 
tection of periodicity (Fig. [8]). Thus, strong consideration 
should be given in favor of including the LMT in a millimeter- 
wavelength VLBI array observing Sgr A* as soon as possible. 
If the parameters of the fully-completed LMT are assumed, 
the bias disappears and the scatter of points on the Hawaii- 
LMT baseline in Figure[5]is similar to that seen on the shorter 
Hawaii-CARMA baseline, and baselines between continen- 
tal North America and the LMT will be of comparably good 
SNR to the lower-resolution PV-PdB baseline. Eventually, the 
LMT and ALMA will be the most sensitive stations in a mil- 
limeter VLBI array and will enable sensitive modelling of the 
Sgr A* system. 

If possible, it would be advantageous to obtain connected- 
element interferometric data of Sgr A* simultaneously with 
VLBI data. While amplitude ratios and phase differences 
on the PV-PdB and (to a lesser extent) SMT-CARMA base- 
lines track the large-scale polarization fraction and EVPA 
fairly well in these models, it is not known what fraction of 
the polarization structure arises from larger-scale emission in 
Sgr A*. If interferometer stations can be configured to pro- 
duce both cross-correlations betweeen telescopes as well as a 
phased output of all telescopes together, opportunities for si- 
multaneous connected-element interferometry may exist with 
the PdB Interferometer, the SMA, CARMA, or ALMA. If sys- 
tem limitations prevent this, it may still be possible to acquire 
very-short-spacing data with those telescopes in CARMA or 
ALMA that are not phased together for VLBI. 

6. CONCLUSIONS 

(Sub)millimeter-wavelength VLBI polarimetry is a very 
valuable diagnostic of emission processes and dynamics near 
the event horizon of Sgr A*. We summarize the findings in 
this paper as follows: 

• Millimeter-wavelength polarimetric VLBI can detect 
changing source structures. Despite low polarization 
fractions seen with connected-element interferometry, 
the much higher angular resolution data provided by 
VLBI will be far less affected by beam depolarization 
and contamination from dust polarization. Polarimet- 
ric VLBI provides an orthogonal way to detect peri- 
odic structural changes as compared with total intensity 
VLBI. 

• Ratios of cross- to parallel-hand visibilities are robust 
baseline-based observables. Short VLBI baselines ap- 
proximately trace the integrated polarization fraction 
and position angle of the inner accretion flow of Sgr A*, 
while longer VLBI baselines resolve smaller structures. 

• Calibration of instrumental polarization terms is not 
necessary to detect a changing source structure, includ- 
ing periodicity, in Sgr A*. 

• Polarimetric VLBI may be able to disentangle the ef- 
fects of rotation measure from intrinsic source polariza- 
tion. Initial results will likely come from observations 
of the timescale of polarimetric variability. If the initial 
array is expanded to allow high-fidelity imaging, polari- 
metric VLBI may be able to map the Faraday rotation 
region and directly infer the density and magnetic field 
structure of the emitting region in Sgr A*. 
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